Magnetic polarons in Ca\- x Y x MnOz 
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Experimental evidence show that in the magnetoresistive manganite Cai- x Y x MnO^, ferromagnetic 
(FM) polarons arises in an antiferromagnetic (AF) background, as a result of the doping with 
Yttrium. This hypothesis is supported in this work by classical Monte Carlo (MC) calculations 
performed on a model where FM Double Exchange (DE) and AF Superexhange (SE) compite. 



PACS numbers: 75.30.Vn, 75.10.-b, 75.40.Mg 

In this work we try to explain experimental resultsB 
on magnetoresistive Ca\- x Y x MnO^ doped manganites, 
that point towards the existence of magnetic polarons. 
The pure sample (x=0) presents G-type AF ordering 
with a transition temperature TV=123K, and insulator 
characteristics. In this sample, all the Mn ions have va- 
lence 4+, where the localized t^g spins have a spin 5=3/2. 
As Yttrium ions are added, itinerants e g electrons with 
spin s=l/2 are introduced and there is a mixed valence 
state of Mn 4+ and Mn 3+ in the proportional — x) and 
x, respectively. The magnetic moment raisesa and the re- 
sistivity drops sharply as x is increased up to x <0.15. A 
thermally activated polaronic behaviour for the conduc- 
tivity was fitted between room temperature and 100K. 
The experimental results were interpreted in terms of 
FM polarons, produced by the alignment of neighbor- 
ing t2 g spins by the e g electrons and immersed in the AF 
background. In this work we compare the experimental 
magnetic measurements with the properties predicted by 
MC calculations on a model where FM-DE and AF-SE 
interactions compete. 
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magnetization, M , was measured with a SQUID magne- 
tometer between 5K and 300K for H <50kOe. 

We measured M vs. H at different temperatures for 
the samples x=0, 0.03, 0.05 and 0.1 and we found a linear 
dependence at high fields: M=Mq+x H. In Fig. 1(a) 
we show Mo vs. T. The behavior resembles that of a 
ferromagnet with a characteristic ordering temperature 
T mo , almost constant for all the samples. For x=0, a 
small M ~0.03^s/f.u. reflects ^a weak ferromagnetic 
Dzialoshynskii-Moriya interaction!! In Fig. 2(a) we plot 
Mq vs. x at T=5K. The curve has an S shape starting 
with a slope of per Mn 3+ ion and reaching a slope 
M Q /x~7fi B at x=0.07. Notice that full FM alignment 
(3^s/Mn site) is not reached. 




FIG. 1. (a) M vs. T, and (b) \ vs - T , for the samples 
x=0 (down triangles), 0.03 (circles), 0.05 (squares), and 0.1 
(up triangles). 

Ceramic polycrystalline samples of Cai-xY x Mn03 
were prepared by solid state reaction methodso. Room 
temperature x-rays diffractograms show that all the sam- 
ples crystallize in an orthorhombic Pnma cell. The dc- 



FIG. 2. (a) Mo and (b) Xi as function of doping x, mea- 
sured at 5K (squares), and calculated (circles). 

In Fig. 1(b), we show x(T) vs. T. The general behavior 
above T mo has been studied in Ref. 1, for £=0, 0.05, 
0.07 and 0.1. The most peculiar feature observed is a 
deviation from the Curie- Weiss law, presenting a negative 
curvature in x _1 (T')- This temperature dependence was 
described in terms of two contributions: 



X{T) = (1 - x)xaf(T) + x \fm(T), 



(1) 



where xaf(T) is the susceptibility of the AF back- 
ground and xfm{T) reflects the strong FM-DE cor- 
relations that results in polaron formation in the PM 
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phase. For x—0, a maximum in %(T) is present at Tjy. 
For x=0.1 the behavior is qualitatively similar: a max- 
imum is observed around T mo , although with a much 
larger value. Below T mo , x(T) decreases rapidly and 
reaches ~ 0.6x(T mo ) at T—5K. For the intermediate 
concentration, x=0.05, a peak at T mo is only a rela- 
tive maximum: xCO increases again below T~70K and 
XAF(5-fT) ~ 1.2x{T mo ). For the low doped x=0.03 sam- 
ple no peak is observed at T mo and x(T) continues to 
increase below T mo up to T~60K, where a broad max- 
imum is present. The large values of x(^mo) take into 
account the progressive polaron formation, indicated by 
the proportionality x(T mo ) oc x, given by Eq.(l), where 
Xaf{T„ 10 ) < XFNi(T mo ). In Fig. 2(b), we show x($K) 
vs. x 

The simplest picture for the behavior of M(H,T) be- 
low T mo , is to associate x(T) with the AF background 
and assume that Mq(T) corresponds to the FM order of 
polarons. However, the measured xCO i s m all cases 
much larger than the CaMnO^ susceptibility. Thus, this 
simple model seems not appropriate to describe the ex- 
periments in the ordered region. In order to obtain a 
better description of the magnetization behavior, we have 
performed a calculation of M vs. H using the numerical 
classical MC technique. 

We consider the following hamiltonian: 

n = - **j (4r c j CT + fc.cj -JhJ2 s *- s * 

{ij)cr i 

+gfx B Hj2S z ,i + Kj2Si-Sj, (2) 

< (»i> 

Here c\ a is the operator creating an itinerant electron 

of spin a at site i, and Sj = X) Q/ 3 c L cr a/3 c i/3 gi ves the 
spin projection of the electron. Due to the strong Hund 
coupling (J n — ► oo) between itinerants e g electrons and 
the localized ti g spins, only parallel spin projections were 
taken into account. In this way, the hopping of itinerant 
electrons depends on the orientation of the localized ti g 
spins, according to the DE expression in one dimension: 
tij=t cos(%/2), where Qij is the relative angle between 
localized spins Si and Sj, t is the hopping parameter. The 
third term, represents Zeeman coupling between mag- 
netic field H and S z ,i, the z component of localized spin 
Si . The AF-SE interaction between localized neighboring 
spins Si and Sj is represented by the last term in Eq.(2). 
Numerical calculations were performed in one dimension, 
using the MC algorithm for a 48 sites chain with open 
boundary conditions. The localized spins were classical 
with modulus 1, and later renormalized to {3+x)fig in 
order to compare with experiments. Local changes in 
t2 g spins were made in conjunction with exact diagonal- 
ization of the itinerant electron system. The resulting 
electronic energy levels were then filled by the available 
number of electrons in the canonical ensemble. We take 
the hopping parameter t=0.1eV as the reference energy, 
and the temperature T=0.005t~6K. Mermin- Wagner's 



theorem precludes magnetic ordering at finite tempera- 
tures in ID systems. However, we estimated a value for 
K from the ordering temperature T mo that would arise 
from a classical Heisenberg model for a system with z=2 
neighbors. For T mo ~100K, we derive KS 2 /t=0.05. 
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FIG. 3. Magnetization M vs. H, for s=0(squares), 
1/48 (circles),.. ,5/48(plus signs), obtained by MC calcula- 
tions. Characteristics error bars are displayed for few points. 

In Fig. 3, we show M vs. H dependence for the 
ID chain, when 0,1, ..,5 electrons are added (equivalent 
to dopings x—0, 0.02, ..,0.1). For the model considered, 
there is a competence between AE-SE and FM-DE inter- 
actions. Preliminary calculationsu show that a polaronic 
phase is obtained for K/t ^0.6 x. These polarons are 
immersed in the otherwise AF background. As we ap- 
ply a magnetic field, the FM polarons tend to align in 
the direction of H . We calculated M(T, x, H) averaging 
over 5000 MC realizations and sites. Mq was obtained 
extrapolating a linear dependence from the highest fields 
(40kOE< H <120kOe). When all polarons are aligned, 
the ratio Mq/x gives a measurement of the polaron size. 
The high field differential slope of M vs. H calculated 
in this case would then arise from two main processes: 
enlargement of the polarons and canting of the AF back- 
ground. 

The slope at high fields is of the same order in all the 
samples. The linear response in the case x—0, is due only 
to the canting of the AF background, x( x = 0) oc K^ 1 . 
For x >0, in the range H <75kOe, we find a rapid growth 
of M vs. H , although strong fluctuations are present due 
to finite-size effects. In this regime the polaron^ origi- 
nally randomly oriented, start to align with HE. The 
numerical results suggests that full alignment is not com- 
pletely reached at 50kOe, which is our experimental limit . 

In Fig. 2(a), we compare the M vs. x dependence ob- 
tained numerically, with the experimental values at 5K. 
As we mentioned above, the slope of Mq vs. x is a mea- 
surement of the magnetic moment of the polaron created 
by each e g electron added. In the experiments, we find 
for x <0.03, that each polaron has a magnetic moment 
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of and for x >0.03, increases up to 7/is/polaron. 

In the present model, we find a linear dependence in the 
curve Mq vs^, i. The size of each polaron is fixed by 
the ratio K/m. The distribution of the itinerant charge, 
shows that each polaron extends over 5 to 6 sites. The 
differential susceptibility x is a measurement of the the 
canting of the AF background and the enlargement of 
polarons with H . In Fig. 2(b) we show the calculated x 
vs. x, that shows an almost constant behaviour. Since 
for :r=0.1, about half of the sites take part in the polaron 
phase, we conclude that the response of polarons to H 
is of the same order to that of the canting of the AF 
background. 

We found a remarkable agreement between experimen- 
tal and numerical values of Mq and x- The numerical 
simulations suggest that linear regime M vs. H, is not 
fully reached at -ff=50kOe, thus experiments at higher 
fields are suggested in order to make a better compari- 
son with numerical prediction, specially in the low doped 
regime, where some discrepancies are observed. 
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